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Abstract 

A simple set of formulas is presented which allows prediction of the fraction of 
ionized carriers in an electron-hole-exciton gas in a photoexcited semiconductor. These 
results are related to recent experiments with excitons in single and double quantum 
wells. 



Many researchers in semiconductor physics talk of "the" Mott transition density 
in a system of excitons and electron-hole plasma, but do not have a clear handle 
on exactly how to predict that density as a function of temperature and material 
parameters in a given system. While numerical studies have been performed for the 
fraction of free carriers as a function of carrier density and temperature lU [2] , these 
do not give a readily-accessible intuition for the transition. In this paper I present a 
simple approach which does not involve heavy numerical methods, but is still fairly 
realistic. The theory is based on two well-known approximations, which are the mass- 
action equation for equilibrium in when different species can form bound states, and 
the static (Debye) screening approximation. In addition, simple approximations are 
used for numerical calculations of the excitonic Rydberg as a function of screening 
length. 

1 Two "Mott" transitions 

The first issue to deal with is what Mott transition we are concerned with. There are 
actually two very different conductor-insulator transitions which go under the name of 
the "Mott" transition. The first, originally envisioned by Mott, occurs when the wave 
function overlap of bound states increases to the point that banding occurs. This Mott 
transition is relevant in the case of doping of semiconductors- when the doping density 
increases beyond a certain point, a mini-band generated from the dopant states will 
arise, in which case the dopant carriers are no longer localized but instead can move 
in Bloch-like delocalized states. The condition for this is ~ a, where a is the 
characteristic length over which the bound state wave function extends, and is the 
average distance between bound states. If we cube both sides of this, we obtain the 
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condition, for a three-dimensional system, 



n~a^ (1) 

where n is the density of dopant electrons. 

This Mott transition condition is equivalent to the condition that the Fermi level of 
the electrons, if they were not in bound states, is much larger than the binding energy 
in the bound states. We can see this by assuming hydrogenic bound states. In this 
case a is given by the Bohr radius, 

where e is the zero-frequency dielectric constant of the medium and m is the effective 
mass of the bound carrier (MKS is used everywhere in this paper). The hydrogenic 
Rydberg for the bound state is 

By comparison, the Fermi energy for a three-dimensional electron gas with two spin 
states, determined by the condition = Jq^^ D{E)dE, where D{E) is the density of 
states, is 

(37r2) ^3 T'/' 
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where n is the electron density. If the Fermi energy is much larger than the Rydberg 
energy, e.g. 10 times greater, we have 
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Essentially, if the Fermi level is large compared to the binding energy in the bound 
states, the binding energy becomes just a perturbation on the Fermi gas. Since the 
Fermi gas is a conductor, the Coulomb binding energy has little effect on that. To put 
it another way, Fermi phase space filling effects, which lead to a Fermi level, become 
important when there is substantial wave function overlap of the electrons, which is 
the same point at which banding occurs. 

By contrast, there is another, thermodnamic insulator-conductor transition which is 
sometimes also called the Mott transition, although Rice has named it the "ionization 
catastrophe." [3] This is usually the relevant transition for photoexcited carriers when 
the excited carrier density is low enough that there is no Fermi level, which is to say, 
the chemical potential is below the ground state energy. In this case, depending on 
the temperature, there can still be a conductor-insulator transition from the carriers 
mostly in excitons to the carriers mostly in plasma. The transition in this case comes 
about due to screening out of the Coulomb interaction which leads to formation of 
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excitons, not due to wave function overlap. The general condition is a ~ where I 
is the screening length, which enters into the screened Coulomb (Yukawa) potential, 
U{r) = g^e-^A/r. The screening length depends on the density; in the low-density 
limit in three dimensions it is given by the Debye screening formula (e.g. 

We can envision this transition in two different directions. If we start with a fully 
ionized, high-density neutral plasma, and reduce the density while keeping the temper- 
ature constant (as, for example, typically happens due to carrier recombination) then 
an instability will occur when the total pair density drops to a point such that the 
screening length is long compared to the excitonic Bohr radius, and exciton formation 
is no longer screened out. Then essentially all the carriers can form into excitons. The 
condition for this is found by reversing the above formula, assuming a is a constant: 

^^^^ (a\ 
n=^^. (6) 

On the other hand, we can imagine starting with a low density gas of only excitons. 
Since excitons are neutral, we do not expect them to contribute to screening, at least 
not until their density is high enough that there is substantial wave function overlap, 
which is the condition ([T| above. However, long before they reach this density, free 
electrons and holes can be generated by thermal ionization. These free carriers can 
leading to screening, which can lead to a reduced exciton binding energy, which in 
turn leads to further ionization. At some critical density this can become a runaway 
process known as the ionization catastrophe, or a thermodynamic Mott transition. To 
determine the condition for this, we must start with the mass-action equation, or Saha 
equation, for the number of free carriers in equilibrium. 



i-ex n>i 



where rig is the density of free electrons (equal to the number of free holes), Uex is 
the density of excitons, and the other factors are the effective density of states factors, 
defined by 
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D,{E)e-^/''^^dE, (8) 



where Di{E) is the density of states for species i, which depends in general on its band 
mass rrii and spin degeneracy. The mass-action equation follows straightforwardly 
from the equilibrium condition that the chemical potentials for the different species in 
a chemical reaction must balance, that is, + = ^^ex^ and from the Maxwellian 
(low-density) formula for the total density, 

n^ = ^-{E.-^^)/kBT^{^^ (9) 

where Ei is the ground state energy of species and from the definition of the binding 
energy Ry = [Ec — E^) — E^x- Number conservation requires that Uex + ne = n, where 
n is the total pair density at any point in time. 



3 



To find the critical density, we set the free carrier density Ue to the value ([6| above, 
instead of the total density. Above this density, screening will make exciton formation 
impossible. Assuming Ry is constant and ng <C riex, so that n ~ Uex, this means the 
critical density for the transition is approximately 



n= , (10) 



where nq = 'i-q^^^q V^q^^- This critical density, for starting at low density and going 
up, clearly gives a very different threshold than ([g]) above, for starting at high density 
and going down. For example, at low temperature, very few excitons will ionize and 
therefore the critical density becomes exponentially large. Of course, at some point the 
Fermi-level Mott transition density ([T]) will be reached. In general, though, the density 
predicted by (10) is well below the density predicted by ([T|. 



The difference between ^ and ( |10[ ) depending on whether the system is starting 
with excitons or neutral plasma leads to the possibility of hysteresis^ that is, the possi- 
bility that the ionized fraction of excitons can depend not only on the total excitation 
density and temperature of the carrier gas, but also on its history. This possibility 
was first presented in Ref. [5] and has been recently revived by Portnoi for two- 
dimensional systems; experimental evidence has been reported in Ref. [7]. 

As seen in this short survey, talking of the Mott transition density is nontrivial. 
The conductor-insulator transition in a carrier gas depends on density, temperature, 
the depth of the bound states, and the history of the system. 



2 Self-consistent calculations for a three-dimensional 
system 

In this section we ignore the possibility of hysteresis and ask for the prediction of 
the mass-action equation for the fraction of ionized carriers as a function of density 
and temperature. If we are to be self-consistent, we must account for the fact that the 
binding energy used in the mass-action equation ([T]) can depend on density through the 
density-dependent screening length. A very useful approximation, based on numerical 
calculations [8j, is 

Ry(/) = |'^='(°>(^-rTlA;)' (U) 

I 0, l/a < 1, 

where I is the screening length. Using this formula, and equation ^ for the screening 
length, we can rewrite ([7]) using number conservation as 

nl = {n- ne)nQe-^^("-^)/^-^, (12) 

which can be solved for rig. 

It turns out that we must adjust this equation slightly. The mass-action equation 
([T]) breaks down in the limit Ry 0, which is exactly the limit we are interested in for 
the ionization catastrophe. This is easily seen by solving for Uex when Ry = 0, which 
gives Hex = nl/nq, which is not zero. We expect, however, that the number of excitons 
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vanishes when the exciton binding energy is zero. The reason for the breakdown is that 
exciton states with extremely shahow binding are stih counted as exciton states, when 
in reahty such states would be unstable, with very large orbits, and electrons and holes 
in such orbits would be indistinguishable in practice from free electrons and holes. As 
a first-order correction, one can subtract from Uex the fraction of excitons which are 
in states with energy greater than Ry, namely a fraction equal to e"^^/'^^^, in which 
case the mass-action equation becomes 

which has the expected behavior that riex — *■ when Ry — > 0. Much theoretical 
effort has been given, especially in the astrophysics community, to justifying this ap- 
proach more rigorously, for example using Planck-Larkin partition functions, which 
give higher-order corrections O |T0] . 

Using this adusted mass-action equation, we then have 

n-' + ni ^ -0 (14) 

If Ry(ne,T) is taken as a constant, then this equation can be solved simply as a 
quadratic equation. This leads to a value for rig which can then be used to recompute 
the value of Ry(ne, T) . This leads to an iterative process to converge on a self-consistent 
value of ng. This iteration converges very rapidly, typically in less than 10 iterations. 
These formulas can be simply programmed in a conditional loop to generate the curves 
shown in this paper. 

Figure 1 shows the fraction of free carriers in bulk GaAs, with the expected behavior 
that the fraction increases toward 100% as temperature is increased and the excitons 
ionize. For simplicity we assume that all holes have the light hole mass; the presence 
of heavy holes will give a correction but not a substantial change to the forms of the 
curves. The dashed curves show the prediction of the mass-action equation for the same 
parameters, but constant excitonic Rydberg. As expected, the self-consistent Rydberg 
solution gives a small correction at low density, and a larger correction at higher density. 
The self-consistent prediction becomes the same as the constant-Rydberg prediction 
at low temperature; for the parameters used in Fig. 1 the two predictions are identical 
below 4 K. 

We can also plot the ionized fraction as a function of density at fixed tempera- 
ture. Figure 2 shows the fraction of ionized carriers as a function of density for various 
temperatures predicted for bulk GaAs. At low temperature, the insulator-conductor 
transition is quite sharp, while at higher temperature it is less so, but at all temper- 
atures there is a "catastrophe" at which rig/n becomes strictly equal to unity. Below 
T = 4 K, the critical density predicted by these equations becomes greater than that 
of the phase-space filling Mott transition density, n~l/a'^ = 3x 10^^ cm^'^. Surpris- 
ingly, the mass- action also predicts that at extremely low densities, the fraction of free 
carriers approaches 100% at all temperatures. This is because the rate of formation 
of excitons is proportional to rig, and at low density the electrons simply cannot find 
each other to form excitons. 

Figure 3 shows the predicted phase boundary. At low temperatures, this is a sharp 
conductor-insulator phase boundary. As seen in Fig. 1, above around 6 K a substantial 
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Figure 1: Solid lines: fraction of free carriers in bulk GaAs predicted by Equation (14) in the 



text, as a function of temperature for the two total pair densities labeled. Parameters used 
are mg = 0.07mo, rrih = O.OSmo, e/eo = 12, and Ry(0) = 4 meV, which implies a = 150 A. 
Dashed lines: the prediction for constant excitonic Rydberg Ry = 4 meV, for the same two 
densities. 
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Figure 2: Fraction of free carriers in bulk GaAs predicted by Equation (14) in the text, as 



a function of density for the four temperatures shown. Parameters for bulk GaAs are the 
same as for Fig. 1. 
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Figure 3: Predicted phase diagram of carriers in bulk GaAs. Heavy line: the critical threshold 
for 100% ionization implied by Equation 14 in the text. Above the boundary, the system 
is a pure plasma. Light line: the phase space filling condition n = l/a'^~3x 10^^ cm"'^. 
Parameters for bulk GaAs are the same as for Figs. 1 and 2. 



fraction of free carriers coexist with excitons even below the ionization catastrophe 
boundary. The effect of the transition wiU be much more dramatic on the luminescence 
spectrum, however, since excitons dominate the luminescence due to their stronger 
oscillator strength. At the ionization catastrophe, the exciton line will disappear. 



3 Calculations for a two-dimensional System 

We can use the same approach for a two-dimensional system, e.g. carriers in quantum 
wells. In this case the screening length is given by ^11] 

1 e^n 

(Note that the Coulomb interaction is still assumed to be three-dimensional, but the 
carriers are constrained to move only in two dimensions.) The two-dimensional exciton 
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Figure 4: Fraction of free carriers in single GaAs quantum wells predicted by Equation 14 
in the text, as a function of temperature for the densities shown. Solid lines: parameters for 
GaAs are the same as for Fig. 1, but Ry = 10 meV a = 120 A, typical values for an exciton 
in a narrow single quantum well. Dashed line: the same parameters, but Ry = 4 meV, for 
the lower density. 

Bohr radius is related to the two-dimensional excitonic Rydberg by tl2j 

In addition, we must use two-dimensional densities of states. A numerical solution 
for the variation of the Rydberg as screening increases in two dimensions has been 
performed by Portnoi and Galbraith [Tl, which can be approximated as 

[ 0, l/a < I, 

where a is the 2D excitonic Bohr radius. 

Figure 4 shows the fraction of free carriers as a function of temperature for the 
parameters of narrow GaAs quantum wells, with Ry = 10 meV, for two densities. In 
this case only one in-plane mass is relevant at low temperature. This figure also shows 
a curve for the same parameters but Ry = 4 meV, with the same exciton Bohr radius, 
which are realistic parameters for indirect excitons in double quantum wells, which 
have been used in recent experiments |13| IT^ . It is perhaps surprising that excitons 
exist all the way up to 100 K or above at these densities. This is a simple consequence 



of the prediction of the mass-action equation (13), which in the hmit of low density 
and high temperature is 

n HQ ksT 
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Figure 5: Fraction of free carriers in single GaAs quantum wells predicted by Equation 14 



in the text, as a function of density for the temperatures shown. Parameters for GaAs are 
the same as for Fig. 4. 



Although higher temperature favors ionization, the free carrier fraction is never 100% 
unless there is an ionization catastrophe. The fraction of the gas which is excitons is 
very low at these high temperatures, (e.g. about 1% at T = 90 K for 4 meV binding 
energy), but since they have a strong oscillator strength they can still have a strong 
emission in the luminescence in many semiconductors. 

Figure 5 shows the fraction of free carriers as a function of total pair density for 
several temperatures, for the case of single quantum wells. A sharp transition occurs at 
all temperatures, above which n^/n is exactly equal to unity. The behavior is similar 
to that of bulk, but with the temperature scale shifted to higher temperature and 
density. Figure 6 shows the predicted phase diagram, for the two different values of 
exciton binding energy discussed above. 



4 Conclusions 

The self-consistent mass action equation presented here allows a relatively simple way 
of calculating the fraction of excitons in equilibrium in a population of excited carriers. 
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Figure 6: Predicted phase diagram of carriers in GaAs quantum wells, for two different 
exciton binding energies. (Heavy solid line: Ry = 10 meV ; dashed line: Ry = 4 meV .) 
Heavy line: the critical threshold for 100% ionization according to the model discussed in 
the text. Above the boundary, the system is a pure plasma. Light line: the phase space 
filling condition n — l/a^. 
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and also the critical density for the thermodynamic Mott transition, or ionization 
catastrophe, in which all excitons are screened out. 

Although the phase boundaries for the thermodynamic Mott transition are not 
simple, their form allows us to speak of two critical thresholds. At moderate densities 
(within an order of magnitude or so of the phase space filling density n = 1/a'^, where 
d is the dimension), one can speak of a critical temperature for the Mott transition, 
since the phase boundary corresponds to roughly constant temperature in this range. 
On the other hand, at high temperature, one can speak of a critical density for the 
Mott transition, since the phase boundary corresponds to nearly constant density in 
this region. This critical density is typically an order of magnitude or more below the 
phase space filling density n = in both two and three dimensions. 

At moderate densities, the critical temperature for the Mott transition is much less 
than Ry/fefi. On the other hand, at low densities, excitons persist as a small fraction 
of the gas to temperatures well above Ry/fcfi. 
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